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Abstract 

We present a simple method of removing the singularities associated with soft photon 
emission to all orders in perturbation theory through exponentiation, while keeping a 
consistent description of hard photon emission. We apply this method to the process 
e + e~ — ► + n7 where we include both Z® and 7 exchange and retain the muon 

mass dependence. The photonic radiation is allowed to be radiated off any charged leg, 
and so we include all initial and final state radiation, as well as all interference effects. 
The effect of exponentiation is to suppress soft photon emission over the cross-section 
you would obtain from working at strictly leading order. We also show how one would 
extend the method to treat the collinear singularity; and remove the associated leading 
mass logarithms. 
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1 Introduction 



Photonic radiation off charged particles plays a very important part in the physics of 
high energy particle colliders. At LEP initial state photonic radiation is responsible for 
shifting the peak and altering its measured width, and hence an accurate description 
of this radiation is necessary to extract the Z boson mass and width in a meaningful way. 
At hadron colliders photonic radiation often forms a background (or indeed a signal) to 
new processes; as, for example, in the future search for an intermediate mass Higgs boson; 
or the testing of anomalous gauge boson couplings. Clearly it is important to understand 
this radiation. 

If we calculate this photonic radiation at tree level in perturbation theory as the 
radiation becomes either soft, or collinear to a massless charged object, we encounter log- 
arithms in E<y/ y/s and respectively. When these logarithms become large the probability 
that additional photons are also radiated becomes large; and so the tree level description 
breaks down. The situation can be improved by going to higher orders in perturbation 
theory, however this does not cure the problem. If we work at next to leading order then 
we include 1 additional photon and so this is an improvement over leading order, but 
when additional photonic radiation becomes important then this next to leading order 
description breaks down, clearly this happens at each finite order in perturbation theory. 
The solution to this is to go to fully infinite order in perturbation theory; and this is 
possible in the soft and collinear limits as in these limits the matrix element can be well 
approximated by the soft approximation and Altarelli-Parisi splitting functions. This 
means that all the logarithms associated with soft emission can be resummed which leads 
to an exponential series [1,2,3,4]. The leading logs or next-to-leading logs associated 
with collinear emission can be resummed and this leads to a "parton distribution" for the 
charged particle [1,5]. These two methods each respectively give an excellent description 
of soft and collinear radiation. The resummation of the soft logarithms can be framed as 
a reordering of the terms of perturbation theory in rigours way as was first done by Yennie 
Frautschi and Suura (YFS) [6]. This method provides an excellent description of both soft 
and hard radiation simultaneously. However there are few Monte Carlo programs that 
incorporate YFS exponentiation and the best (YFS3[3],BHLUMI4.0[4]) only describe 1 
hard photon exactly (with a 2 nd included through the leading log splitting function). On 
the other hand tree level Monte Carlos at a fixed order perturbation theory can describe 

- 1 - 



arbitrarily many hard photons [7, 8] . In this paper we take our previous tree level Monte 
Carlo [8] that calculated e + e~ — > + rry and exponentiate the soft photons. We use 

a less rigorous, but more simple, form of exponentiation than full YFS exponentiation. 
Our form of exponentiation is equivalent to YFS exponentiation except that the effects 
of the virtual loop Feynman diagrams are not explicitly included, but only added using 
an ad hoc method. We show how one would go about resumming the collinear radiation 
to remove the large mass logarithms that occur, however as our primary interest in this 
paper is the interface between hard and soft radiation we make no attempt to include 
this resumed collinear radiation in our Monte Carlo program. Now when we calculate 
radiation at tree level we are forced to include cuts that keep us separate from both the 
soft and collinear region so we do not encounter the singularities that are present there; 
in fact, as we have suggested above these cuts should be large enough that we do not 
approach the areas of phase space where the soft and collinear logarithms become large. 
Now in our work as we resum the soft logarithms and remove the soft singularity, and 
so we are no longer required to keep the cut that keeps us separate from the soft corner 
of phase space, however we are forced to retain the cut that keeps us separate from the 
collinear singularity as we have not dealt with the singularities that are there. As with 
tree level calculations this cut should be viewed as an experimental cut, that should be 
chosen large enough to keep us away from the large collinear logarithms, that is imposed 
on all photons. 



2 Exponentiation of soft photons with other hard photons present 

In any process in which charged particles are accelerated the soft photon approxima- 
tion tells us that soft photons are predominately radiated off external legs, and that the 
matrix element for the Feynman diagram with a soft photon radiated off external leg % 
with charge e is, 

M e (2.1) 
k-pi 

where M. is the matrix element for the process without a photon, and e and k are respec- 
tively the photon polarisation vector and momentum. 

So if the matrix element for the process, 

-> e+(pi) e"(p2) V + (P3) H~{pa) + rvy (2.2) 
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is M. n the matrix element squared for the process, 

-> e+(pi) e _ (p 2 ) Ai + (p 3 ) A*~(P4) + ™7 + 7*(*0 (2-3) 
in the soft photon approximation is given by, 

l^-WK-^t^-^) 1 (2.4) 

\ pi • k p2-k P3 • k p^-k J 

after we sum over the different spin states of the photon, where we definite all particles 
in the final state. This means that initial state particles will have negative energy. 

If the soft photon and the n original photons are in mutually exclusive areas of phase 
space then there is no symmetry factor between the two. We can force this to happen by 
only considering the original n photons with energy, E^, larger than some cut, and soft 
photons with energy, E s , less than that cut, i.e., 

E h > E cnt > E s (2.5) 

then differential cross-section is given by, 



2 1 i KA ,2 ( PI P2 , P3 P4 



2 



da n+1 = - e z — \M n \ z ^ - ^ + ^ - ^L- rf(LIPS) n+3 (2.6) 
Flux \Pl'k P2-k ps-k p^-k J 

\Mn\ 2 d{UPS) n+2 (Eikonal Factor) (2.7) 



where, 



Flux 1 1 v '""^ v ' 2E(2n 



(Eikonal Factor) = -e z ( ^- + ^- ) (2.8) 

P\-k P2-k p$-k P4-k 1 



'2 



/(«) 



(2.9) 



and integrating over the soft photon phase space for photon energies E m [ n < Ej < E cu t 
gives, 

(E cut \ 



1 , . . |2 at 



d<r n ,l = ^ | Mnr ^(LIPS)n+2 ^ ln (^7^ J ^(^) (2-10) 



n 



with, 



«?(ft c ) = f f(n)<m (2.ii) 



2 We give the explicit form of g in the appendix. 
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and so we see from (2.10) that the probability to emit a soft photon just factorises the 
lower order cross-section. However in the factorising term we have a contribution from 
ln(E m [ n ) and this diverges as E m [ n goes to zero, this is the soft singularity. 3 

Now for fixed particle momenta in da n unitarity tells us that when we integrate over 
the soft photon momenta that da n \ is finite. As the soft, non colinear, singularity contains 
no electron or muon mass terms this tells us that the soft logarithm must cancel in the full 
calculation including virtual diagrams[9]. We can achieve the same effect as the virtual 
diagrams by imposing an effective lower energy cut off on the photon energy, -Ereg, and 
for the higher order corrections corrections to be 0(a) we require, 

In = 0(1) (2.12) 

where \fsj1 is the maximum photon energy. So the cross- sect ion, differential in n photons 
with energy larger than E cu ^ and integrated over 1 soft photon with energy smaller than 
E CVLt is given by, 

In doing this we improve the accuracy of the differential cross-section calculation from 
(D(ahi(E)) to 0(a). The lack of knowledge of E reg beyond (2.12) is exactly the lack 
of knowledge that we have from doing a tree level calculation without calculating the 
virtual diagrams. Also note that if we consider the differential cross-section for values of 
E cn t < E re g then this appears to go negative ! This is just an artifact of working, at least 
as far as the soft logs go, beyond tree level. 

If we now consider additional soft photonic radiation then for the process, 

0^e+(pi) e~(p 2 ) li + (p 3 ) H~(p4)+ni + Js(h)---ls(km) (2.14) 

then the soft photon approximation gives the matrix element squared as, 

m 



i=l 



3 Also note that the term g(Q c ) diverges in the limit of massless fermions when the photon 
becomes collinear to any fermion. This is the collinear singularity, isolating this singularity 
and removing it leads to collinear description of photon radiation. In this work we will 
apply angular cuts to all photons (be they hard or soft) to keep the photons separate 
from the fermions, and as such we work in an area of phase space where these collinear 
singularities are not important. 



The phase space for these soft photons also approximately factorises, 

m E^dE' dVt ■ 

rf(LIPS) 2+n+TO ~ rf(LIPS) 2+n J] 2 ^ ( ^ )3 ? (2.16) 

and so as before we can write the differential cross-section for the process with n photons 
with energy larger than E cu t and integrated over m soft photons with energy less that 
£cut as, 



rn 



da n , m = da n - In ^— J g(n c ) j (2.17) 

where the 1/m! term is the symmetry factor that we get from integrating the m identical 
soft photons over the same regions of phase space. 

If we now ask what the differnetial cross-section, rfcr^ xp , for process where we produce 
an arbitary number of photons satisfying the cut Q c , where we stay differential in n 
photons with energy larger than E cu t and integrate over an arbitary number of photons 
with energy less than E cu t, we find, 

do"n XP = da n fi + do n i + da n ^ 2 + ■ ■ ■ 

= <*» f 1 + (£ (ft ) H + h(& ln (ft) *H 2 + ■ 



(1*) 

V -^reg / 
/ EcuA M^c)/4^ 2 ) 



da„ exp ( ^ln (Jgj g(Q c , 




\ -^reg / 

<i(T^ xp has the sum over an infinite number of soft photons, and this means that it has 
all the real soft logarithms resummed in it, each with the soft singularity removed; this 
means it will provide an accurate description of arbitrarily soft photons. Notice that in 
this "exponentiated" form we have no problems with cross- sect ions going negative for 
values of E cv x < E TC g. da n * p has no knowledge of photons that fail the angular cut Vt c 
and with respect to those photons <ia^ xp is a strictly leading order quanity. This means 
da^ does not describe this collinear radiation, and in particular when large collinear 
logs appear in g(£l c ) or da n , da n xp will not give a good description of the radiation. 

Now we can use the form of (2.21) to arrive at a workable Monte Carlo stratergy; 

• Choose the number of photons n to be produced. 
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• Calculate a point in phase space for the process (2.2) where the n photons have energy 
larger than E cn \,. 

• Calculate the hard matrix element for the process (2.2) and weight it by, 



• Sum over the events to perform the integration over n and the hard phase space for 
the process (2.2) . 

Notice that we can only calculate the expression (2.22) after calculating the point in 
the hard phase space as the g{VL c ) term depends upon the orientations of the charged 
particles in the process. 

In this prescription we appear to have introduced two extra parameters, E re g and 
E cnt , over the tree level calculation. Now the value of E Te g is introduced to cancel 
the soft singularity; and as in a tree level calculation we don't know how much of the 
singularity cancels, the value of E re g is unknown beyond what we learn from equation 
(2.12). Thus the dependence of our calculation on E Te g is a measure of our uncertainty in 
only calculating the real photon emission diagrams without calculating the virtual photon 
loop diagrams that cancel the soft singularities. The same level of uncertaintly is also 
present at tree level although there is no parameter that displays it. 

Our apparent dependence upon E cu ^ is more worrying. _E cu t was just introduced as 
a parameter to distinguish the "hard" photons from the "soft" photons, however there is 
no physical meaning to the words "hard" and "soft"; there is no magical energy where 
photons suddenly become "hard". Consider what happens as we decrease E cn i; for the 
hard matrix elements we approach the soft singularity and the cross-section grows rapidly; 
now this increase in cross-section comes about from low energy photons, as these photons 
have low energy they do not change the physical signature of the process. However as the 
hard matrix elements are growing the eikonal factor (2.22) decreases and if the value of 
E cu t is to have no physical significance then these two effects should cancel each other. 
This will give us a very strong test of our results, if we have a result independent of E cn t 
we are practically required to a correct calculation of both the hard matrix element and 
the soft eikonal factor. Notice that as we include an arbitrary number of soft photons in 
the soft exponentiating factor we are required to consider an arbitrary number of hard 
photons in order for the E cn t cancellation to occur. This is true even if those hard photons 




(2.22) 
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would not be experimentally observed on energy grounds alone. This is the opposite of 
what we do in a strictly LO calculation; there if we require say the 2 photon differential 
cross-section it is essential that we generate exactly 2 photons, here it is essential that 
we also calculate the n photon rate for all n > 2. Also notice that we will only get this 
cancelation between the soft photon exponentiation and the hard photon cross-section if 
we treat both the soft and hard photons identically, in particular this means that the 
angular cuts that we apply to the hard photons must also be applied to the soft photons. 

Comparing the naive form of exponentiation considered here with the more rigorous 
YFS exponentiation, then YFS exponentiation explictly includes both hard photon cor- 
rections and virtual photon corrections to the exponentiating factor. The virtual photon 
corrections are totally ignored in this work, as we consider no virtual diagrams; this means 
that this work has uncertainties of 0(a) ~ 1% - these uncertainties manifest themselves 
as a lack of knowledge in the value of E re g. The hard photon corrections give only a finite 
contribution to the matrix element in the soft corner of phase space, and as the volume 
of this phase space tends to zero as E cu t — > in this limit the hard photon corrections 
have vanishing effect. In practice this means that we should take -E C ut small enough that 
our results have no dependance upon -E C ut- in YFS exponentiation the total energy of all 
radiated soft photons is forced to be less than some cut; whereas in our form of exponen- 
tiation it is each soft photon energy individually that has energy less than some cut; this 
means that there is a worry that although each soft photon individually only carries away 
a small amount of energy, as there are in infinite number of soft photons they may carry 
away a sizable amount of energy in total. If this were to happen then our results would 
be incorrect because energy would not be conserved in the E cu t — > limit; however if this 
were to happen our results would not be independent of the cut E cu t, thus independance 
of our results upon this cut check that our energy cut on the soft photons is valid. 

3 Collinear radiation 

In this paper we are primarily interested in the interaction between soft and hard 
radiation, as such we have only included the effects of resumming soft radiation in our 
Monte Carlo program. This means that we are always forced to include an "experimental" 
angular cut on all photons, Q c , that keeps us separate from the singularity when radiation 
becomes collinear to a massless charged object. As we have this "experimental cut" our 
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method has nothing to say about photons that fail the cut Q c , or the physics that these 
photons generate (like the shift in the measured Z° mass and width). In this section we 
show how such collinear radiation could be included in our method at the leading collinear 
log level; however we make no attempt to include this collinear radiation in a Monte Carlo 
set up. 

In principle we would like to proceed in a similar way to the soft radiation, that is 
we evaluate the matrix element for a charged particle radiating an arbitrary number of 
collinear photons, then integrate the photon momenta over the collinear region of phase 
space; and finally sum over an arbitrary number of collinear photons. In practice the 
matrix element for a charged particle radiating an arbitrary number of collinear photons 
is not known and so this method fails. This means that we can not calculate all real 
collinear logarithms in the way that we calculated all real soft logarithms; however we 
can calculate the leading logs or the next to leading logs through an evolution equation. 
If we have 2 collinear photons then if one photon is far more collinear than the other we 
can approximate the matrix element as two independent collinear emissions of photons. 
Usually the evolution equation is in terms of the virtuality of the initial of final charged 
particle, or the maximum that the radiation can have; however for our case it is far 
more convenient for the evolution to be done in terms of the angle of the emitted radiation. 

If we have an massless charged particle, p, that radiates a collinear photon, k, at an 
angle 9; so the final charged particle, p' , has a fraction z of the initial particle energy, 
then the lowest order matrix element is multiplied by, 



where, 



and in the collinear region the phase space is given by, 

d(LIPS) = ^ dz d{2p'.k) = ixE'E 1 dz dcos6 (3.3) 

where we have integrated over the unimportant azimuthal angle. So the differential cross- 
section to emit an extra single collinear photon is given by, 



\M spUt \ 2 d(LIPS) = (27re 2 ) d{- ln(l - cos &)) V(z)dz (3.4) 
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If we now consider the cross-section to emit an arbitrary number of collinear photons up 
to some angle 6 then the lowest order cross-section is multiplied by a function, D; for this 
to be useful we need to know the energy fraction of the final charged particle, z, and so 
it is convenient to define, 

V(z) = ^ where D = [ dz V(z) (3.5) 
dz J 

then the evolution of T>{z) is given by, 

dV(z) 



= ( 27 re 2 ) / —V(x)V(z/x) (3.6) 
Jz x 



d(Zn(l -cos0)) 

Where to derive this equation we have assumed that the splitting function for several 
angular ordered photons is given by a production of splitting functions, this is strictly 
only true when the photons are strongly ordered. This means that this evolution equation 
only sums the leading collinear logarithms. 

If we integrate (3.6) over z we find, 

D = A(l -cos#) 2 ^ e2p (3.7) 

where P = Jq dz V(z). Usually within the + prescription P = 0, and so D is constant 
when evolved in 9. This means that as we vary the angular cut arround the charged 
particles Q c the cross-section in the collinear region does not change. However the cross- 
section for the noncollinear radiation does depend upon the angular cut Q c , and this 
means that the total cross-section is not independant of this cut. In order to restore 
indepenence of the cross-section on fl c we define V , 

rl-E leg /E 

T(z)=V(z)+5(z-l) dyV(y) (3.8) 



where E is the energy of the charged particle, and then use V in place of V . 

Equation (3.6) being an evolution equation only tells us how radiation up to some 
angle from the charged particle is related to radiation at any other angle, it does not tell 
us where to start the evolution. Now if the charged particle were indeed massless then the 
starting point of the evolution would be uncalulable, as in QCD parton distributions for 
the proton and fragmentation functions; however for massive charged particles radiation 
is supressed in a dead cone surrounding the charged particle of angle defined by sin 6>rjc — 
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m/E, where m and E are respectivly the mass and original energy of the charged particle. 
This suggests the starting point of the evolution as, 



V{z,& DC ) = 5{z-\) (3.9) 

This gives that, 

D(e BC ) = 1 (3.10) 

r2 N 2ire 2 P 



f9F 2 \ 

i.e., A ~ (^t) , this means that, 



,2\ 2^e 2 P 



2E Z 

D(cos6 « 0) ~ I — s- ) (3.11) 



m 2 



Now in the limit m — > 0, D diverges; this, at least for final state radiation, is unphysical[9] . 
The origin of this logarithmic divergence is in the integral over cos 6? in the differential 
cross-section, this takes the form, 



5 o 



cos9^y/l-m 2 /E 2 /2F 2 \ 

d(-ln(l-cos0) ~ln (— o-l (3.12) 
cos0wO \m l J 

This is the collinear logarithm that, like the soft logarithm, also cancels (for final state 
radiation [10]) on the virtual diagrams. As with the soft singularity we have not calculated 
the virtual diagrams we again do not know how much of the mass logarithm cancels, and 
so we again introduce an extra parameter that quantifies our lack of knowledge in this 
cancellation. For the mass logarithms to be absent in the total cross-section we require, 

D(cosQ w 0) = 1 + 0(ct) (3.13) 

and for this to be true we choose A in equation (3.7) to have the value, 

X — l + fia where /i ~ 1 (3-14) 

This means that, 

, 2 \ 27r e 2 P 

Dfo c ) ~ (1 + /ia) ^ j (3.15) 
and so we choose the starting point for the evolution (3.6) to be, 

/ o \ 2 ne 2 P 
I 171 \ 

V(z, 9 BC ) = S(z - 1) (1 + (^) (3.16) 



This ensures that the total cross-section is not enhanced by logarithms in m. For angles 
less that 9-qq the fragmentation function V no longer evolves in cos 9 but stays at the value 
of V{9y>q). It should be noted that we require a different value for \i for initial and final 
state radiation; indeed whereas for final state radiation we know that /i ~ 1 [10] , we know 
no such thing for initial state radiation where the cross-section may contain logarithms in 
the mass. Also note that T>(z) plays a very different role for initial state radiation than 
for final state radiation. For initial state radiation the charged particle fragments into 
a charged particle and photons before the hard scattering, and so the charged particle's 
momentum is degraded and the hard scattering takes place at a lower y/s\ whereas for 
final state radiation the fragmentation of the charged particle happens after the hard 
scattering and so does not change the hard scattering. 

Having set up the fragmentation functions D and T> the hard differential cross-section, 
where we have summed over an arbitary number of soft and/or collinear photons, is given 
by 



, col— exp , 
dcr n * =da n 



Zl,Z 2 

/Ecut\M^V 47r2 ) (3.17) 
V -^reg / 

x V(Q C , z\) dz\ V(Q C , z<i) dz2 T>{VL C , z%) dz% V(£1 C} z±) dz^ 



x 



where da r 



is the hard differential cross-section where the incomming electrons have 
fractions z\ and z^ of the beam energy; z% znd z^ are the fractional energy that the final 
state muons have. 

As the soft radiation is radiated off particles with a similar virtuality as the collinear 
radiation, or alternativtly they are radiated at the same timescale, it is not clear wether 
g(f2 c ) should be calculated before of after the collinear radiation has been emitted in the 
T> functions. However recall that g is the integral over the matrix element squared, where 
the matrix element is the sum of terms that go like p/ p-k; and this depends only on the 
direction of p and not its energy in the massless limit. So g only depends on the directions 
of the charged particles, and not on the amount of energy they carry. This means that 
g is independent, in the massless collinear limit, of wether it is calculated before or after 
collinear radiation is emitted. 

With the differential cross-section in the form (3.17) it is easy to give a Monte Carlo 
mechanism for generating hard photons. 



• Choose the number of photons n to be produced. 

• Choose z\ and z 2 and hence the initial electron energies. 

• Generate the hard phase space for n photons where the initial electrons have fraction 
z\ and Z2 of the beam energy. 

• Choose z% and z\ that give the fractional energy that the observed muons have over 
the muons that emerge from the hard scattering. 

• Calculate the hard matrix element for hard phase space and weight it by, 



• Sum over events to perform the integration over the hard phase space, z\ through Z4 
and n in the usual Monte Carlo way. 

In the soft case the dependance upon E cn t cancels exactly between the hard cross- 
section and the soft exponential factor. In this collinear case we do not expect exact 
cancelation of Q c dependance between the collinear factor D and the hard cross-section, 
as we have only resummed the leading logs in D rather than all the logs as in the soft 
case. Non the less there should only be rather mild dependance on Q c arrising from the 
sub leading logs that occur in the hard cross-section, but not in the factor D. 

The resummation of these collinear logs has no analogy in YFS exponentiation, which 
makes no attempt to resum collinear logs. Within YFS exponentiation the full collinear 
logs can only be added as hard corrections to the exponentiating factor. This means that 
our collinear resummation should give a more accurate description of radiation in the 
collinear region than YFS exponentiation. 

4 Numerical Results 

We have applied the procedure from section 2 to the process, 




v(n c , zi ) v(n c ,z 2 ) v(n c ,z 3 ) v(n c ,z 4 ) 




for arbitrary n. To do this we need the hard matrix element for (4.1), this we have 
calculated in a previous paper [8] . In that calculation we retained the full mass dependence 



of the muons, and allow an arbitrary number of hard photons radiated from both the initial 
and final state fermions. We use the approximation that the electron is massless, however 
as we have angular cuts to keep the photons separated from the electrons we do not work 
in an area of phase space where the electron mass terms are important. We include both 
s channel Z° and photon exchange. 

For the physical constants we use, 

4m = V128 
«em = 1/137 

sin% = 0.23 (4.2) 
M z o = 91.175 GeV 
m M = 105.6584 MeV 

where we use a em (M z ) = a4m as the coupling for the s channel and photon, and 
a em (0) = af^f as the coupling for all external photons [11]. 

As we have not removed the large logs associated with collinear radiation in the Monte 
Carlo we are required to keep away from the regions of phase space where these logs be- 
come important, this means that we must keep all photons (be they exponentiated or not) 
separate from all charged particles. This we do by imposing the following experimental 
cuts, 

> 5° 

(4.3) 

| cos# 7 | < 0.9 

where 9^ is the angle between a muon and photon, and 9^ is the angle of a photon from 
the beam pipe. In order for the final state particles to be observed we impose, 

9^ > 20° 

35° < cos0 M < 145° (4.4) 

where 9^ is the angle between the two muons and 9^ is the angle of a muon from the 
beam pipe. We choose, 

= M z (4.5) 

Now for large n the hard matrix element for process (4.1) can be very time consuming to 
calculate, so we calculate up to 4 photons exactly - then for larger numbers of photons 
we calculate the 4 most energetic photons exactly, and then use the soft approximation to 
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Fig.l The contributions from varying numbers of hard photons to the 
inclusive cross-section for the process e + e~ — > as a function of E cn t. 

We choose E ie g = 10 GeV, \fs = M z and use the other cuts described in 
the text. On the right we expand the scale and just show the inclusive 
cross-section. We also show the Monte Carlo errors associated with each point. 

calculate the remaining photons. In practice the contribution to the cross-section from 4 
or more hard photons is always very small, and usually several of the photons are very soft; 
this means that we make only a small error by using this approximation, while speeding 
up the code considerably. 

The first test of the program is the dependence upon E CVLt and E Te g. In Fig.l we 
show the contributions from the n hard photon cross-sections to the total inclusive cross- 
section as a function of -E C ut- Clearly these each individually have a large dependence 
upon -Ecut- However when we sum the contributions from each number of photons to 
form the physical total inclusive cross-section where photon radiation satisfies the cut fl c 
we obtain the solid line which clearly has far less dependence upon E cn t. We also show 
the total inclusive cross-section with an expanded scale. It is clear that for small i?cut the 
cross-section is independent of E cn t within the Monte Carlo error; however there is also 
a clear rise in the total inclusive cross-section for values of E cu ^ ~ 1 GeV. This can be 
understood as the soft approximation that we have used to derive the exponential factor 
(2.22) breaks down for photon energies larger than Oil GeV). 

As we take E cv & smaller it takes a longer time to calculate the cross-section with 
the same accuracy as for smaller E cv ^ values we have a larger contribution from larger 
numbers of hard photons; and the cross-sections for large numbers of photons are very 
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E re g(GeV) 

Fig.2 The dependence of the inclusive e + e~ — > exponentiated 
cross-section on E Te g. We also show the LO result for the same cross-section. 

time consuming to calculate due to the rapid increase in the number of Feynman diagrams. 
We have also calculated the E cn t dependence of the more exclusive 1 and 2 photon cross- 
sections, and for small E cn t values there is no visible dependence upon E cu t. This gives 
us confidence that our computer code has no mistakes in it. In all the results that follows 
we will choose E cu ^ = 10~ 2 GeV which is safely in the region where our results are 
independent of E cu t, also the value of E cu t is not so small that we spend large amounts 
of time calculating the cross-section for large numbers of hard photons. 

If Fig.2 we show the dependence of the inclusive cross-section upon E Te g, and also the 
LO e + e~ — > cross-section. The dependence of the cross-section on E rC g is clear 

to see. The two cross-sections are not directly comparable as the exponentiated result 
has dependance upon Q c ; however it is clear that the cross-sections are comparable as we 
expect, and the difference is of 0(a) for values of E ve g ~ 10 GeV as suggested by equation 
(2.12). In all following results we will choose E re g = 10 GeV. 

Moving onto the more exclusive 1 photon differential cross-section, if we consider 
events with 1 or more photons present then the energy of the most energetic photon is a 
physical quantity. This means that we can form the differential cross-section da/dlnEi, 
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Fig. 3 The 1 photon differential cross-section 
da /din E\ plotted both at LO and exponentiated. 

where E\ is the energy of the most energetic photon, in a meaningful way. If we do this 
we find the results shown in Fig.3; in this figure we have also plotted the LO prediction 
for da/d\nE\. As E — > we see that the da L0 /d\nE\ — > constant, this is manifest from 
the soft photon approximation. Equation (2.10) tells us that, 



da\° 
din Ei 



a 



(To — ~g(Q c ) ~ constant 



(4.6) 



That this breaks down for E\ ^ 1 GeV we can understand as being due to initial state 
radiation (ISR), the large cross-section comes from the process e + e~ — > Z — > with 
an extra photon radiated; however if the photon is radiated from the initial state with an 
energy greater that 0(Y %) this forces the exchanged boson far off mass shell, and this 
suppresses the cross-section. 

We can see that for large E\ the LO 1 photon cross-section gives an answer similar 
to the exponentiated 1 photon cross-section; however for smaller E\ the LO result over 
estimates the cross-section. The effect of resumming the logarithms associated with the 
soft singularity is to change the shape of the 1 photon cross-section. Notice that changing 
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FigA The ratio of the exponentiated divided by the LO 1 photon 
differential cross-section as a function of the photon energy, E\. 

the value of E re g largely just changes the overall normalisation by a factor of, 

(£ r e g r a5( " c)/4 ^ 2 

and does not change the overall shape of the 1 photon cross-section. 

For E\ ~ \fsj1 the cross-section has a small maximum, this is mainly due to e + e~ — > 
production followed by a muon fragmenting into a photon. It is clear to see why 
this produces a maxima for large photon energies, consider, 



(4.7) 



M2) 7(3) 



(4.8) 



this fragmentation is proportional to the propagator of muon 1; i.e., it is inversely pro- 
portional to p\ = (p2 + P3) 2 , the final muon photon invariant mass. If we look at a 
fragmentation at a particular p\ = (p2 +P3) 2 value, if the photon and muon 2 have sim- 
ilar energies then they are produced largely collinear and the cut (4.3) vetos this event. 
However if either the final muon, 2, or the photon is very soft then the angle 6^ is far 
larger and we pass the collinear cut (4.3). The configuration where the photon carries 
most of the energy of muon 1 gives rise to the small peak in Fig.?, at large E\. 
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In Fig A we plot the ratio, 

da L0 /d\nEi V ' ; 

We can see that there is no particular energy where exponentiation becomes important, 
rather thats as we consider lower and lower energy photons exponentiation slowly becomes 
a more important factor. For our choice of -Ereg = 10 GeV the exponentiated result is 
comparable to the LO results for large E\ (the dip close to y/s/2 is due to the the process 
(4.8) mentioned before), however for photon energies of 1 GeV the exponentiated result 
is down to 90% of the LO result. 

It is of interest to ask how the soft photon approximation compares to the exact matrix 
element result. However before we compare the two approaches there is a subtlety with the 
soft photon approximation that should be appreciated. The soft photon approximation 
tells us that the matrix element with a photon present is the matrix element without that 
photon multiplied by an eikonal factor. However what does the matrix element without 
the photon mean? Consider the process e + e~ — > /z + /i~7, if we remove the photon then 
we no longer have 4 momenta conserved; in particular if we ask what is the Q 2 carried by 
the s channel Z° or 7 we can form this in two ways as (p e + + p e -) 2 or as (p^+ + £>^-) 2 
and these have different values. As such the soft photon approximation has ambiguities 
associated with what the matrix element without photons is. Clearly as the photon energy 
goes to zero this ambiguity disappears. 

In the current case of e + e _ — > with y/s = M% this ambiguity can be of great 

importance for photon energies greater than 0(Tz)- If we have initial state radiation 
(ISR) with energy greater than 0(Yz) then this forces the exchanged boson off mass 
shell, and this strongly suppresses the cross-section; this means that the relevant Q 
for the exchanged Z® is (p^+ + £> )U -) 2 which means that for photon energies greater than 
0(Yz) the Z® propagator is strongly suppressed. Similarly if we have final state radiation 
(FSR) with energy greater than OiT %) then the relevant Q 2 is (p e + + p e -) 2 - With this 
in mind, and given that photon radiation is radiated off a pair of charged legs (at least 
for massless fermions) a good prescription is to boost the two legs that the photon can be 
radiated off in such a way that 4 momentum is conserved and that their momentum only 
changes by a small amount. In the present case this is messy as the soft eikonal factor 
(2.4) has 6 different combinations of legs that it can be radiated off, and so there are 6 
different boosts that need to be carried out. This is beyond the scope of the current work; 
instead we consider the two simpler cases of pure ISR and pure FSR (i.e., the cases where 
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Fig. 5 The soft photon approximation compared to the hard matrix 
element, both exponentiated, for the cases of pure ISR and pure FSR. 

we take the final and initial legs respectively to be chargeless). For pure ISR (FSR) it is 
clear which legs the photon is radiated off, and hence which legs should be boosted; for 
ISR we boost the electron legs, and for FSR we boost the muon legs. If we are to boost 
momenta p\ and P2 we use the boost defined by, 



/ 

Pi 



ail PI + anP2+PT 
021P1 + &22P2+PT 



(4.10) 



with p'y = pf, p'2 = p\ PvPt = an d P2'Pt = 0j where we require p'^ and p^ to satisfy 
4 momentum in the system without photons. 

We plot the ratio of the soft approximation to the hard matrix element in Fig.5. It 
is clear that the soft approximation is good for photon energies up to O(10 GeV) but 
underestimates the hard result for larger energies. For the case of ISR this makes little 
difference as radiation with energy larger than 0(T^) is strongly suppressed; and so we 
only make mistakes where the contribution to the cross-section is small. The rate for hard 
FSR can be quite appreciable and for this hard radiation the soft approximation underes- 
timates the cross-section by a large factor. It seems hard to improve the implementation 
of the soft photon approximation in such a way that it correctly describes these hard 
photons. 

As a final numerical result we wish to consider the photon photon invariant mass 
distribution in the process, 

+ /i + /i- 77 (4.11) 
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Fig.Q The exponentiated differential cross-section da/dM^ both with hard 
matrix elements and in the case where the hardest photon is treated exactly 
and the remaining photons are described by the soft photon approximation. 

this process has been the subject of much interest since the L3 collaboration at LEP 
announced 4 events with M 77 = 60 GeV [12] in apparent conflict with the Standard 
Model. Although the conflict with the Standard Model has since gone away with improved 
statistics this process serves as a good illustration of a process that contains both soft 
and hard photons and so both exponentiation of soft photons and exact matrix elements 
for hard photons are important. 

We use cuts similar to the L3 experimental cuts; 

|cos# 7 | < 0.9 36° <9p< 144° 

# 7M > 5° 6m > 20° (4.12) 

E 1 > 1 GeV 

We show the exponentiated hard differential cross-section we obtain in Fig.Q. Now we 
can consider various different approximations to this process. If we just use the LO 
hard matrix element result we would expect to describe large M 77 well, but have larger 
difficulty with small M 77 ; in practice with the cuts used there is no difference within the 
Monte Carlo errors from the exponentiated hard result for all values of My 7 . 
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We also compare the soft photon approximation, for the comparison we have boosted 
the final state muons in order to conserve 4 momenta within the soft approximation. This 
keeps the exchanged on mass shell even when we have hard ISR, as such we would 
expect this to overestimate the cross-section. We show the exponentiated cross-section 
where we calculate the hardest photon using the exact matrix element and calculate 
additional photons using the soft approximation. It is clear that although this works 
well for small M 77 this underestimates the cross-section by a large amount for large 
M 77 . Surprisingly the differential cross-section is almost unchanged when we calculate 
all photons within the soft approximation from the result when we treat one photon 
exactly. 

For this process it is clear that describing the hard photons correctly is more important 
than exponentiating the soft photons. 

5 Conclusions 

In this paper we have presented a method for exponentiating the singularities associ- 
ated with soft photon emission while still retaining exact matrix elements for hard photon 
emission. We have applied this method to the process, 

e + e~ — > ji + jjT + n'-f (5.1) 

where the photons can be radiated from any of the charged legs, i.e., we have included 
ISR, FSR and all the interference terms. We have also retained the mass dependence of 
the muon. 

Exponentiation is of importance for low energy photons where it suppresses the cross- 
section over the LO prediction, whereas exact matrix elements are important for hard 
photons when the soft approximation that is the basis of exponentiation breaks down. 
Experimentally an accurate description of both soft and hard photons is often simulta- 
neously required as in a process where the hard photons are of primary interest the low 
energy photons are used to normalise the theoretical predictions for that process. 

We have give a method by which one could also remove the collinear singularities, 
and resum the associated logarithms, in section 3. As we are not currently interested in 
these singularities we have made no attempt to include these effects in our Monte Carlo; 
but instead impose experimental angular cuts on all photons. 
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Appendix 



In this appendix we give the explicit form of the function g(Q c )- 9 is defined as, 

g (n c ) = [ f(n)dn (a.i) 

with, 

Now /(Jl) contains two types of term, P 2 p\-p2 I Pl'k P2'k and E 2 m 2 / {p\-k ) 2 . Concen- 
trating on the first of these we have to integrate, 



h = j 2E 2 ^ yz , <m 



I (Ei 



2 PVP2 
pi-k p2-k 

2(E 1 E 2 -Pi P2 cos p) d cos 9 



(A.3) 



P\ cos 9) (E2 — P2 cos p cos 6 — P2 sin p sin 6* cos 0) 
where we have written, 

Pl = (Pi, Pi, 0,0) 

P2 = (E 2 , P 2 cos p, P 2 sin p, 0) (A.4) 
k = (E, E cos 9, E sin # cos 0, P sin sin 0) 

Integrating over gives, 

f 27r dcos9 , , . 

j., — I ( A 5 1 

7 (Ei - Picos9)((P 2 cos9 - E 2 cosp) 2 + m 2 sin 2 p) 1 / 2 
and performing the indefinite integral over cos#, 

h(pi,P2,COS0) = —— 2 ^ 2 d2 m/2 X 



In 



(y 2 + m 2 sin 2 p) l l 2 (z 2 + m 2 sin 2 pP 2 ) 1 / 2 + zy + m| sin 2 p Pi 



Pl - Pi cos 9 

(A.6) 

where y and z are defined as, 

y = —P 2 cos 9 — E2 cos p 

(A.7) 

z = P2E1 — E 2 pi cos p 

Now we want integrate f(Q) over the region defined by fi c , that is the whole of space 
with small angular regions removed about the charged particles. Two of those charged 
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particles are p\ and p2, now 9 is the angle between p\ and the photon, and so we can 
form the integration over the region with > 9\ and #27 > #2 as ? 

I = (h(Pi,P2,cos6 1 ) - h(pi,P2,-l)) - (h(p2,Pi,l) - h(P2,Pl,cos6 2 )) (A.8) 

Now as m2 — > 0, Ii(pi,p 2 , — 1) and Ii(p2,Pl, 1) diverge, and these two divergences cancel 
on each other. In order to avoid numerical difficulties it is wise to extract the singular 
terms and cancel them by hand. This we can do by rewriting the logarithm in equation 



and cancel them by hand. 

So far we have assumed that the angular cuts about charged particles do not overlap, 
in practice if the angular cuts overlap the charged particles are close in phase space and 
so p\-p 2 is small and the eikonal factor as a whole is small. Numerically the problem 
that we encounter is that I goes negative, in the case that this happens we set I = and 
as the eikonal factor is small this only makes a small error. 

In addition to particles 1 and 2 we also have two additional particles 3 and 4 that the 
photon is kept separated from. However these particles are typically far away from the 
singularities of / and is slowly varying. This means we can approximate the effect of the 
cuts #3 7 > #3 and 9^ > 64 by evaluating E 2 p\-p<i / p\-k p 2 -k in the direction of particles 
3 and 4 and then multiplying this by the volume of space that the cut excludes, 



(A.6) as, 





and 




we then subtract these two contributions from (A.8). 



It remains to calculate the integral of E 2 m /(pi-k) . This is given by 



/JT^fYl^ r r2ir ^2 
-77 d£l = l d cos 9 / d6- -77 
(Pi-k) 2 J Jq V (Ei-Picos9) 2 



(A.ll) 



2?rm 2 



Pi (£1 -Pi cos 9) 



As before we impose the angular cut about particle 1 by integrating (A. 11) between 
— 1 < cos 9 < cos#i, and impose the angular cuts about particles 2,3 and 4 by evaluating 
E 2 m 2 ' /(pi-k) 2 in the directions of particles 2,3 and 4 and multiplying by the volume of 
the angular cut. 
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